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Abstract 

The consistency of the concept of quantum (quasi)particles possessing effec- 
tive mass which is both position- and excitation-dependent is analyzed via 
simplified models. It is shown that the system may be stable even when the 
effective mass m = m(x, E) itself acquires negative values in a limited range 
of coordinates x and energies E. 



1 Introduction 



Non-relativistic quantum dynamics of point particles is most often studied 
via ordinary differential Schrodinger equation 
ft 2 d 2 

— y(x) + V(x)V(x) = Ey(x), *(i)6L 2 (K) (1) 

2m ax 2 

where the real function V(x) characterizes an external local potential while 
the particle mass m > is just a given constant. Various generalized forms 
of this equation were introduced due to the practical needs of the description 
of motion of a particle or quasiparticle inside a medium. The medium may 
make the mass position-dependent [TJ. The phenomenological appeal of the 
spatial variability of the effective mass is accompanied by some additional 
formal merits of the generalization, say, in the context of supersymmetric 
quantum mechanics [2]. 

In the majority of papers which studied the models with m = m(x) the 
authors assumed that m{x) > [3J. The more we were surprised when 
we found [I] that in the so called VT— symmetric version of the quantum 
Coulomb problem the stability of the system required an anomalous, asymp- 
totically negative choice of the coordinate-dependent effective mass. This 
observation re- attracted our attention to quantum models with m = m{x) 
and, in particular, to their subset in which one encounters the anomalous 
m(x) < 0, in some nonempty interval of coordinates at least. 

In what follows we intend to describe some results of our analysis. Firstly, 
in section [2] we introduce an elementary solvable Schrodinger equation with 
a piecewise constant mass such that m = m(x) < for x G (—1,1). We shall 
demonstrate that the bound-state spectrum of energies of such a model is 
unbounded from below so that the model cannot be interpreted as realistic 
due to its immanent instability. This result supports our expectations that 
for the models with indefinite mass m{x) such an instability will be generic. 

One of the possible methods of elimination of similar pathologies is then 
proposed in section [3j We recall the general Feshbach's definition of the ef- 
fective quantities in quantum theory (cf. ref. [5]) and conjecture that the 
apparent anomalies in the behavior of the benchmark square-well model of 



2 



section [2] are not realistic. We argue that their emergence should be at- 
tributed to our unfounded complete suppression of the necessary variability 
of the effective mass with the energy. On this basis we propose and de- 
scribe our first exactly solvable model with m = m(x, E) in section [3] and its 
alternative version in section HI 

Our main result is that both of these amended benchmark models may 
remain stable (and, hence, acceptable for phenomenological purposes), pro- 
vided only that the mass m(x,E) stays merely anomalous (i.e., negative) in 
a restricted range of its arguments (i.e., in finite intervals of coordinates x 
and energies E). 

A few formal aspects and possible consequences of the reinstalled math- 
ematical consistency and acceptability of Schrodinger equations with sign- 
changing effective masses m(x, E) will be finally mentioned in section El We 
shall re-emphasize there the phenomenological as well as purely theoreti- 
cal appeal of the use of effective masses m(x, E) which change their sign. 
Whenever it happens just locally in both x and E, one may encounter no 
contradiction with expectations and/or with the general principles of quan- 
tum mechanics. 

2 The instability of quantum systems with a 
locally negative effective mass m = m(x) 

In the majority of applications one accepts the requirement of the positivity of 
the mass m(x) > as natural. In the light of Ref. [I], such a rule might have 
its exceptions. We should admit that in the latter paper the weird-looking 
asymptotic negativity of the effective mass can in fact be attributed to our 
rather technical postulate of the loss of the observability of the coordinate 
x in the asymptotic region [6j. Vice versa, we believe that the effective 
mass m(x) should remain asymptotically positive whenever the asymptotic 
coordinates x remain real. In this context, our present letter may be read as 
motivated by the question of possible existence of some physically consistent 
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scenarios admitting a negative effective mass m(x) < 0, say, in some not too 
large finite interval of coordinate i£l 

For the sake of definiteness let us simplify the mathematics of such a 
conceptual problem by considering just the toy model in which the energy 
spectrum is discrete and in which one works with the deep square-well po- 
tential 

, i +00, \x\ > L > 1 , 
V(x) ' (2) 

10, \x\< L. 

Moreover, in our models the mass will be piecewise constant, say, as follows, 

*.(,) = { l ^ (1 ' L) - (3) 

I ™>o, \x\ <1 . 

Unfortunately, whenever mo is negative, the spectrum of the bound-state en- 
ergies of such a model becomes unbounded from below. This makes the whole 
system unstable with respect to perturbations so that its physical meaning 
becomes highly questionable. Still, it makes good sense to understand this 
counterintuitive and apparently discouraging fact in a more quantitative de- 
tail. 

Let us consider the toy-model Schrodinger equation 

with potential (jSJ) and mass (j3J) where, say, m = — 1. Let us also restrict 
our attention, for the sake of brevity, to the mere even-parity bound-state 
solutions %f(x) = $f(—x) such that $f(—L) = *(L) = 0. 

In units h 2 /2 = 1 and in the first step of analysis we shall assume that 
E = k 2 is non-negative. This means that we have to solve two differential 
equations, 



-^^>(x) = k 2 ^(x) xe (-L,-l), 

-£^(x) = -k 2 ^(x) xe(-i,o), 



with the respective explicit solutions 

#(x) = sin k(x + L) x e (-L, -1 



(6) 

= Dcoshkx ie(-l,0) 



(?) 



which satisfy the usual physical boundary conditions. In addition, these 
solutions must also satisfy the usual mathematical matching conditions at 
x = —1, 

D cosh k = sin k(L — 1) , 
— D sinh k = cos k(L — 1) . 
This determines the eigenvalues E n = k\ via an elementary transcendental 
equation 

tanh/c x tanfc(L — 1) = — 1 (8) 
(cf. the two samples of its graphical solution in Fig. [1]). 




Figure 1: Graphical solution of Eq. (jSJ) at L = 2 (the small circles mark the 
intersections of — tanh/c with cot k(L — 1)). 

In the same units h 2 /2 = 1 and in the second step let us set E = —k 2 
and write down the two corresponding equations, 



x e (-L, -l) 
a; 6 (-1,0), 



—£i*(x) = k 2 V{x) 
with the respective explicit solutions 

ty(x) = sinh k(x + L) 
ty(x) = F cos KX 

The matching conditions at x = — 1 yield the two constraints, 

sinh k(L — 1) = F cos k , 
cosh k(L — 1) = F sin k . 



(9) 



x e (-L,-l), 
x e (-1,0) . 



(10) 



ID 



The resulting implicit definition of the negative-energy eigenvalues E n = —k^ 
is then obtained, 

tan k x tanh k(L — 1) = 1 (12) 

and yields infinitely many roots again (cf. their graphical representation in 
Fig. [2]). Hence, the spectrum appears unbounded from below. This forces us 
to declare the model manifestly unphysical. 




Figure 2: Graphical solution of Eq. ( II 2p at L = 2 (the small circles mark the 
intersections of tank with cothfc(L — 1)). 

The key role of our particular example ([3]) possessing the discrete bound- 
state spectrum which is not bounded below should be seen in its methodi- 
cal importance. Indeed, although the model is mathematically correct, its 
stability will be immediately lost under the influence of virtually arbitrary 
perturbation. 

Naturally, such a feature of our toy model is generic. Moreover, its insta- 
bility under perturbations might have been expected a priori since, after the 
kinetic energy becomes negative, the perturbed system will certainly have a 
tendency of plunging, spontaneously, into the lower and lower energy states. 
The corresponding wave functions will acquire a highly oscillatory form. 

Qualitatively, the perturbative instability phenomenon will remain the 
same for virtually any external potential V(x). The exact solvability of our 
model gives such a sensitivity to perturbations just a fully explicit form. 
For our particular constant and energy-independent choice of the parameter 
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mo < one observes the explicit form of the steady increase of the localization 
of the unperturbed wave functions inside the non-empty interval of x. With 
the growth of parameter k we merely witness the emergence of an explicit 
square-well-like part of the bound-state spectrum which is turned upside 
down. 



3 A realistic, energy-dependent model 

In the language of physics, the unperturbed versions of our simple-minded 
m < toy model should be re-interpreted as incomplete. Any realistic 
realization of similar systems (which, in principle, radiate, i.e., act as a source 
of energy) must necessarily be reconsidered as coupled to an environment. 
In other words, we must re-classify our physical Hilbert space as a mere 
subspace "Hp of a larger physical Hilbert space r Hj vXl . 

The lack of the knowledge of the (usually, prohibitively complicated) full 
Hamiltonian H* which would characterize the environment (and which 
would act in the full Hilbert space T-L^ ul1 ) leads to the necessity of acceptance 
of some hypotheses. Fortunately, the use of the formalism of the Feshbach's 
projection-operator techniques (cf. [5] for details) appears extremely efficient 
in this context. 

For our present purposes, in particular, is it sufficient to take into account 
that even if the rigorously known spectrum of our unperturbed, decoupled 
(or, in the Feshbach's language, P— projected) Hamiltonian of section __ is 
not bounded from below, any of its realistic perturbed versions may be given 
the Feshbach's semi-explicit form 

H (effective) = p jjfull p + p jjfull q _ ______ Q ft full p ^3) 

where Q = I — P. Such an effective Hamiltonian may certainly represent 
measurable phenomena but, by assumption, our knowledge of this operator 
is restricted and incomplete. Its important merit is that it is still defined in 
the accessible, "small", P— projected Hilbert space Tip. Next, the energy E 
need not be considered complex in the present setting. The reason is that 
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we may and shall tacitly assume that the whole spectrum of H* remains 
discrete and that, for the sake of simplicity, the energy E in ( 1131) does not 
belong to the spectrum of Q H* Q. 

The only conclusion which we can make about the difference between 
the schematic, unperturbed negative-mass Hamiltonian (say, P H^ ul1 P) and 
the whole family of its possible realistic effective descendants ( fl3l) is that 
all of the latter operators must be, necessarily, manifestly energy-dependent. 
This fact may be recalled as giving the reasons why the choice of the effective 
mass should also be considered energy-dependent in general. Vice versa, once 
we admit that the mass is energy dependent, m = m(x,E), any restriction 
of attention just to the positive values of this phenomenological parameter 
becomes entirely artificial. 

Via our previous, energy-independent example we already understand 
that the emergence of the unlimited decrease of the sequence of the bound- 
state energy levels should be attributed to the purely mathematical role 
played by the negative constant m < at the large parameters E = —k 2 . 
On this basis we may most simply stabilize the spectrum by keeping the 
effective mass positive beyond certain threshold, m(x,E) > for E < E t hr- 

In parallel, the physical meaning of the threshold energy cut-off may vary 
with the hypotheses concerning the environment. The choice of its value 
may be interpreted as a compressed information about the interplay between 
the complicated coupling operators P H^ ul1 Q and the Hamiltonian of the 
environment Q H' Q. Naturally, even such a weak form of information 
about the hidden and/or prohibitively complicated dynamical mechanisms 
may still play the role of a phenomenological source of variability of the 
effective mass. 

Needless to add, the use of the energy-dependent effective operators finds 
a broad range of applications in various domains of quantum physics [Tj . The 
key feature of their Feshbach's mathematical origin is that every effective op- 
erator (defined in the P— projected subspace) varies with the changes of the 
energy of the (P + Q)— projected system. In this sense, the unwanted emer- 
gence of mathematical anomalies (like, typically, the unbounded spectrum as 
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mentioned above) may still be eliminated and attributed to an inappropriate 
treatment of the energy-dependent simulation of the effects of the environ- 
ment. 

For the sake of definiteness of our argument let us now proceed by teach- 
ing by example again, replacing the constant parameter m in Eq. fl3]) by 
its suitable energy- dependent generalization. It will be allowed negative at 
E > E t hr for, say, E t hr = 0. For the sake of simplicity, let us also add the 
convenient assumption that mo(E) — > 1 for E — > — oo and that m (E) — > — 1 
for E — > +oo. 

Under these assumptions, we may select, for illustration purposes, the 
following, most elementary interpolation ansatz replacing Eq. ([3]), 

{1, \x\ E (1, L) , 

II V , J, ^ 
-tanh (E), \x\ < 1 . 

Under this choice the value of the mass parameter just simulates the emer- 
gence of an anomaly in the kinetic-energy operator for a restricted range of 
the energies. As already mentioned, such an effective kinetic-energy oper- 
ator finds its natural hypothetical origin in the Feshbach's reduction of a 
realistic or "complete" (i.e., formally, (P + Q)— projected) Hilbert space (in- 
cluding some unspecified and formally eliminated "medium") to its suitable, 
explicitly tractable (i.e., P— projected alias model-space) subspace. 
In the resulting amended version 

--;±-=r-fz9(x) = E*(x), *(-L) = ¥(L) = (15) 
m{x, E) ax z 

of our toy-model Schrodinger equation the emergence of energy-dependence 
of the mass certainly does not change the applicability of the matching 
method of solution. The insertion of the amended effective mass (|14p will 
still enable us to proceed in an almost complete parallel with the preceding 
section. First of all we shall split our equation in its "outer" and "inner" 
part, 



-£ I *(x) = E*(x) xe(-L,-l), 
^y(x) = -Etanh(E)V(x) xe(-l,0) 
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and restrict our attention, for the sake of brevity, just to the even-party states 
again. Next we distinguish between the positive and negative sign of E. In 
the former case we set E = k 2 yielding the ansatz 

^(x) = sink(x + L) , x e (-L, -1) , 



$>(x) = DcoshXx , are (-1,0), A = \(k) = k Vtanh k 2 > 0, 

(17) 

while in the latter case we put E = — k 2 and get 

\&(x) = sinh k(x + L) , x G (— L, —1) , 



\I/(x) = F cosh/xa; , x G (— 1,0), /i = /i(k) = k\/ tanh /t 2 > . 

(18) 

The respective matching conditions read 

D cosh A(/c) = sin k(L — 1) , 

w v ; (19) 

-X(k) D sinh A(A;) = k cos fc(L - 1) , 

sinliK(L — 1) = Fcoshu(K) , 

v ) " ; ' (20) 

k cosh k(L — 1) = — /i(/t) F sinh //(/«) . 
They imply that the eigenvalues may be obtained from the respective tran- 
scendental equations 

Vtanh k 2 tanh \(k) tan k(L - 1) = -1 , (21) 



vt^mh/t 2 tanh h(k) tanh k(L — 1) = — 1 . (22) 

We see that the left-hand side of the latter equation is nonnegative for any real 
k so that the set of its roots is empty. The bound-state energies themselves 
are all determined by Eq. f l2T|) (cf. Fig. [3] which may be perceived as just a 
slightly deformed analogue of Fig. [TJ. 

The even-state part of the spectrum proves bounded from below. This 
is our main conclusion. Along similar lines, the same conclusion may be 
obtained for the odd-parity bound states. We may summarize that whenever 
one wants to work with the quantum systems in which the effective mass 
can get negative, one is not allowed to neglect the variability of the effective 
mass with the energy. 
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Figure 3: Graphical solution of Eq. (f2Tj) at L = 2. 

4 Another exactly solvable model 

The introduction of the energy-dependence in the effective mass has been 
based on the existence of a hypothetical environment. One would expect, on 
this basis, that the "realistic" ^—dependence of m(x, E) would be smooth. In 
such a case, the study of the related ordinary differential effective Schrodinger 
equations may proceed in parallel with the E— independent cases. A compact 
review of these parallels may be found in Ref. [8]. It has been emphasized 
there that the overlaps between the wave functions corresponding to differ- 
ent energies will not vanish in general. In our present particular model the 
readers might easily check this fact using the available explicit formulae for 
wave functions. 

An exhaustive clarification of this apparent paradox may be found in 
Ref. [8]. Interested readers find there not only the standard Hilbert-space 
interpretation of the energy-dependent Hamiltonians (based on the use of 
suitably adapted ad hoc inner products) but also the explicit realization of 
the bi-orthogonality and bi-orthonormality relations between eigenstates in 
similar models. One of the most concise reviews of the extension of this 
formalism to the case of the general non-Hermitian observables may be also 
found in Ref. [9]. 

We saw, via the schematic illustrative example of preceding section, that 
the bound-state spectrum might be very sensitive to the value of the threshold 
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energy Ethr- In particular, it may be expected much more sensitive to this 
value than to the concrete shape of the function w,q{E) itself. For this reason, 
let us now replace the smooth-function ansatz (1141) by the simplified but 
arbitrarily shifted step-function shape 



m(x, E) 



\x\ G (1,L) 

-1, E>E thr I , , (23) 

+1, 




where the threshold energy is now to be chosen as any finite, real and, for the 
sake of definiteness, negative variable E t h r — ~P 2 - 111 this context one must 
be aware of the fact that for a generic phenomenological potential given in 
advance, the similar discontinuities and abrupt changes of the mass might 
give rise to singular contributions to the effective potential. This would 
deserve a separate analysis. Interested readers are recommended to check 
the existing literature in this respect [HI]. 

Certainly, there will be no solutions at the energies E < —(3 2 . For the 
sake of brevity, we shall again discuss just the even-parity bound states. In 
this case, the half-infinite interval of the eligible energies should again be 
split into the high-energy half-line with E = k 2 > and the complementary 
low-energy interval with E = — h 2 where one must merely admit the limited 
range of the eligible k G (0, (3). 

In the former case we may return to Fig. Q], reminding the readers that the 
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high-energy subset of our present, "second-example" eigenvalues E n = k, 
n = 1, 2, . . . will still be defined as roots of the same, unchanged Eq. (JHJ). 

The similar implicit specification of the second, low-energy subset of the 
remaining eigenvalues E\^ n = k^, n — 1, 2, . . . , N should again be deduced, 
mutatis mutandis, from the unmodified secular Eq. (fT2l . The problem is 
easily resolved since a return to Fig. [2] reveals that we just have to replace 
the original infinite half-axis of k of section [2] by its finite subinterval of 
admissible k G (0, (3). 

We may conclude that the new lowest (i.e., ground-state) energy level 
will always emerge when there appears a new root of Eq. ([12"]) . i.e., the new 
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root k n = fi CT iticai of equation 

tan /3 critica i x t&nh [3 critica i(L - 1) = 1 . (24) 

In other words, with the increase of the real variable (3 > 0, new and new 
elements of the sequence of the energies E\- n = —^n computed in section [2] 
and sampled in Fig. |2] will be reclassified as the acceptable elements of the 
spectrum of the present new model. 

Special attention must be paid to the limiting case in which k = Perineal- 
It is necessary to keep in mind that also in such an extreme case the matching 
condition is satisfied in standard manner so that the related solution has to 
be accepted as a valid new lowest-lying bound state. Perhaps, it is worth 
adding that as long as, by our assumption, the effective mass of such a state 
remains negative in the whole interval of x G (—1,1), one should not be 
surprised that such a state is more or less fully localized in this interval. We 
also see from the explicit formula fflUj) for wave functions that the number 
of the nodal zeros of such a wave function, paradoxically, increases with the 
decrease of E thr = -f3l itical . 

5 Discussion 

Naturally, our present, methodically motivated restriction of the number of 
the free parameters to the necessary minimum might be very easily relaxed 
in any future work. One might point out, for example, that the level-pattern 
as provided by Fig. [2] is, via secular Eq. ffT2]) . closely related to the role 
and interpretation of the threshold-energy parameter Ethr — ~P 2 entering 
our second illustrative example of paragraph HI This means that the same 
qualitative picture will be also provided by the trivially rescaled models while 
the pattern may be changed by an introduction of an additional parameter. 

In such a setting it has been extremely interesting for us to notice that the 
replacement of our one-parametric toy model (GJ) + ([3]) by its two-parametric 
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generalization with optional a > and b = b(E) > in the effective mass 

m(x) J \ (25) 
-1/6 , \x\ < a 

leads to the mere replacement of our previous secular Eq. ([12]) by the rescaled 
relation 

hhCL 

tan — x tanh k(L — a) = b . (26) 

Naturally, no qualitative changes are encountered for the single rescaling 
variable b = a. Nevertheless, a nontrivial qualitative change of the pattern 
will emerge in the "deep narrow mass well" regime with small a and b such 
that also the ratio a/b := v gets small. 

In this regime, secular Eq. f)26p does not possess any small— k solution so 
that tanh k(L — a) = 0{1). One reveals that in the way illustrated by Fig. |2l 
the sequence of the "larger" roots k 2 ~ irb/a, k 3 w lixbja "disappears" to the 
right infinity in the limit v = a/b — > 0. The remaining, single and numerically 
easily obtainable leftmost root K\ will remain finite. Its approximate value 
will be given by the reduced secular equation 

K\ = — coth K\L . (27) 
v 

Obviously, the right-hand-side expression does not change too quickly with 
K\ and it will be, in addition, not too much larger than the ratio bjv itself. 
Thus, typically, for bjv = 1 one obtains, numerically, the approximate value 
K\ ~ 1.2 of the root. 

One just has to add that in such a special b = v limit our mass function 
(1251) with —1/b 2 ~ 1/a may be interpreted as the Dirac's delta function so 
that the survival of just the single bound state in the spectrum might have 
been, intuitively, expected. 

The most important mathematical merit of our present choice of the 
position-dependent effective masses (which are just piecewise constant func- 
tions of x) is that the selection of the rigorously defined self-adjoint Hamilto- 
nian is traditional and trivial. In the first two single-parameter toy-models it 
was provided by the most elementary matching of the logarithmic derivatives 
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of the wave functions. This observation has been complemented by the third 
example in which we revealed that one could also employ some less trivial 
matching prescriptions in certain limiting cases. 

Still, no use of the full-fledged self-adjoint-extension theory [TT] was needed. 
Due to the piecewise constant nature of our effective masses, we did not 
even have to consider their usual von Roos' p2] factorizations m(x, E) = 
mi(i, E) rri2(x, E) 777.3(2;, E), nor did we have to select an appropriate order- 
ing of the individual mass factors and differential operators d/dx (which, 
naturally, do not mutually commute in general). 

The resulting simplification of the mathematics and of the necessary func- 
tional analysis did certainly make the physical interpretation of our models 
more transparent. This enabled us to emphasize that the Feshbach- method 
origin of the concept of the effective mass opens several new perspectives in 
the related phenomenology. In our present letter we managed to demonstrate, 
first of all, that, and under which conditions, the effective mass m = m(x, E) 
need not necessarily be required a positive function of x and/or E. 

We might point out that our present detailed description of a few toy 
models clarified that there exists an intuitively obvious connection between 
the violation of the positivity requirement m(x) > and the breakdown of 
certain traditional theorems. Let us recall, for example, our observation that 
in the cases of indefinite m(x), the ground-state wave function is allowed to 
possess nodal zeros so that the traditional Sturm-Liouville oscillation theo- 
rems cease to be valid and must be modified. In parallel, we noticed that in 
many cases, the anomalous increase of the number of the nodal zeros with 
the decrease of the bound state energy is accompanied by the perceivable 
localization of the wave function inside the interval of negativity of m(x, E). 

All of these mathematical observations might find their potential future 
physical applications in all of the phenomenological models where one has to 
mimic some effects of an (unknown, implicitly described) external medium 
by means of the use of the (energy-dependent) effective operators. In this 
context we kept in mind some possible parallels with the relativistic phenom- 
ena and with the well known contrast between the behavior of particles and 
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antiparticles. Thus, we concentrated our attention to the toy-model study of 
the extreme scenario in which the effective kinetic energy reverses its sign. 

We may conclude that the oversimplified nature of our illustrative exam- 
ples served our purposes well. In particular, they helped us to clarify the 
connection between the decrease of the value of the negativity-threshold en- 
ergy E t hr and the increase of the number of the "anomalous" , localized and 
quickly-oscillatory low-lying bound states. Via our second example admitting 
a freely variable E t h T — —ft 2 -, this connection has been even quantitatively 
sampled by the elementary-looking restriction of the pattern given by Fig. [2] 
to the finite interval of admissible k < j3. 
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